Running methratio on files that have been trimmed for adapters

In [5]:
#working directory (parent)
wd="/Volumes/web/Mollusk/bs_larvae_exp"

#where is bsmap
bsmap="/Users/Shared/Apps/bsmap-2.73/"

#Where is bsmap file
bsmapfile="/Volumes/web/Mollusk/bs_larvae_exp/AdaptorTrimmed_BSMAP_files/"

#Output for methratio file
methratio="/Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_methratio_files/"

#genome file 
genome="/Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa"
In [6]:
!python {bsmap}methratio.py -d {genome} -u -z -g -o {methratio}methratio_out_M1.txt -s {bsmap}samtools {bsmapfile}M1/condor-stdout-0.sam
@ Thu Apr  3 13:39:42 2014: reading reference /Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa ...
@ Thu Apr  3 13:40:15 2014: reading /Volumes/web/Mollusk/bs_larvae_exp/AdaptorTrimmed_BSMAP_files/M1/condor-stdout-0.sam ...
[samopen] SAM header is present: 7658 sequences.
	@ Thu Apr  3 13:45:56 2014: read 10000000 lines
@ Thu Apr  3 13:46:26 2014: combining CpG methylation from both strands ...
@ Thu Apr  3 13:46:55 2014: writing /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_methratio_files/methratio_out_M1.txt ...
@ Thu Apr  3 13:57:22 2014: done.
total 8716599 valid mappings, 48672250 covered cytosines, average coverage: 1.78 fold.

In [7]:
!python {bsmap}methratio.py -d {genome} -u -z -g -o {methratio}methratio_out_M3.txt -s {bsmap}samtools {bsmapfile}M3/condor-stdout-0.sam
@ Thu Apr  3 13:57:23 2014: reading reference /Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa ...
@ Thu Apr  3 13:57:56 2014: reading /Volumes/web/Mollusk/bs_larvae_exp/AdaptorTrimmed_BSMAP_files/M3/condor-stdout-0.sam ...
[samopen] SAM header is present: 7658 sequences.
	@ Thu Apr  3 14:04:07 2014: read 10000000 lines
@ Thu Apr  3 14:05:08 2014: combining CpG methylation from both strands ...
@ Thu Apr  3 14:05:37 2014: writing /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_methratio_files/methratio_out_M3.txt ...
@ Thu Apr  3 14:18:44 2014: done.
total 9773350 valid mappings, 53390397 covered cytosines, average coverage: 1.78 fold.

In [8]:
!python {bsmap}methratio.py -d {genome} -u -z -g -o {methratio}methratio_out_T1D3.txt -s {bsmap}samtools {bsmapfile}T1D3/condor-stdout-0.sam
@ Thu Apr  3 14:18:46 2014: reading reference /Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa ...
@ Thu Apr  3 14:19:19 2014: reading /Volumes/web/Mollusk/bs_larvae_exp/AdaptorTrimmed_BSMAP_files/T1D3/condor-stdout-0.sam ...
[samopen] SAM header is present: 7658 sequences.
@ Thu Apr  3 14:22:28 2014: combining CpG methylation from both strands ...
@ Thu Apr  3 14:22:57 2014: writing /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_methratio_files/methratio_out_T1D3.txt ...
@ Thu Apr  3 14:29:21 2014: done.
total 5759274 valid mappings, 26507549 covered cytosines, average coverage: 1.32 fold.

In [9]:
!python {bsmap}methratio.py -d {genome} -u -z -g -o {methratio}methratio_out_T1D5.txt -s {bsmap}samtools {bsmapfile}T1D5/condor-stdout-0.sam
@ Thu Apr  3 14:29:22 2014: reading reference /Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa ...
@ Thu Apr  3 14:29:56 2014: reading /Volumes/web/Mollusk/bs_larvae_exp/AdaptorTrimmed_BSMAP_files/T1D5/condor-stdout-0.sam ...
[samopen] SAM header is present: 7658 sequences.
@ Thu Apr  3 14:34:49 2014: combining CpG methylation from both strands ...
@ Thu Apr  3 14:35:18 2014: writing /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_methratio_files/methratio_out_T1D5.txt ...
@ Thu Apr  3 14:45:50 2014: done.
total 6974302 valid mappings, 45446985 covered cytosines, average coverage: 1.54 fold.

In [10]:
!python {bsmap}methratio.py -d {genome} -u -z -g -o {methratio}methratio_out_T3D3.txt -s {bsmap}samtools {bsmapfile}T3D3/condor-stdout-0.sam
@ Thu Apr  3 14:45:52 2014: reading reference /Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa ...
@ Thu Apr  3 14:46:47 2014: reading /Volumes/web/Mollusk/bs_larvae_exp/AdaptorTrimmed_BSMAP_files/T3D3/condor-stdout-0.sam ...
[samopen] SAM header is present: 7658 sequences.
	@ Thu Apr  3 14:57:00 2014: read 10000000 lines
@ Thu Apr  3 14:57:38 2014: combining CpG methylation from both strands ...
@ Thu Apr  3 14:58:31 2014: writing /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_methratio_files/methratio_out_T3D3.txt ...
@ Thu Apr  3 15:19:06 2014: done.
total 8848034 valid mappings, 52256465 covered cytosines, average coverage: 1.65 fold.

In [11]:
!python {bsmap}methratio.py -d {genome} -u -z -g -o {methratio}methratio_out_T3D5.txt -s {bsmap}samtools {bsmapfile}T3D5/condor-stdout-0.sam
@ Thu Apr  3 15:19:08 2014: reading reference /Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa ...
@ Thu Apr  3 15:20:02 2014: reading /Volumes/web/Mollusk/bs_larvae_exp/AdaptorTrimmed_BSMAP_files/T3D5/condor-stdout-0.sam ...
[samopen] SAM header is present: 7658 sequences.
	@ Thu Apr  3 15:30:41 2014: read 10000000 lines
@ Thu Apr  3 15:31:16 2014: combining CpG methylation from both strands ...
@ Thu Apr  3 15:32:09 2014: writing /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_methratio_files/methratio_out_T3D5.txt ...
@ Thu Apr  3 15:50:44 2014: done.
total 8808529 valid mappings, 51732717 covered cytosines, average coverage: 1.69 fold.

In [1]:
#NOTEBOOK NAME
nb='Methratio_AdapterTrimmed_Larvae'

#NOTEBOOK LOCATION
nbloc='/Users/claireolson/Dropbox/Lab\ \(1\)/IPython_nbs/Claire/'



#TIMESTAMP
date=!date +%m%d_%H%M
#LOCATION OF POST DIRECTORY
post='/Volumes/web/Mollusk/iPython_Posts/'
In [2]:
cd {post}
[Errno 2] No such file or directory: '/Volumes/web/Mollusk/iPython_Posts/'
/Users/claireolson/Dropbox/Lab (1)/IPython_nbs/Claire

In [3]:
!ipython nbconvert {nbloc}{nb} --to html --template full --output {nb}{date}
[NbConvertApp] Using existing profile dir: u'/Users/claireolson/.ipython/profile_default'
[NbConvertApp] Converting notebook /Users/claireolson/Dropbox/Lab (1)/IPython_nbs/Claire/Methratio_AdapterTrimmed_Larvae.ipynb to html
[NbConvertApp] Support files will be in Methratio_AdapterTrimmed_Larvae[0404_0946]_files/
[NbConvertApp] Loaded template html_full.tpl
[NbConvertApp] Writing 207630 bytes to Methratio_AdapterTrimmed_Larvae[0404_0946].html